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Charm production in association with an electroweak gauge boson at the LHC 
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The production of charm quark jets in association with electroweak gauge bosons at the LHC 
can be used as a tool to constrain quark parton distribution functions (PDFs). Motivated by recent 
measurements at the Tevatron and LHC, we calculate cross sections for W/Z + c, comparing these 
to W/Z + jet, for various PDF sets. The cross-section differences can be understood in terms of 
the different underlying PDFs, with the strange quark distribution being particularly important 
for W + c production. We suggest measurements of appropriately defined ratios and comment on 
how these measurements at the LHC can be used to extract information on the strange and charm 
content of the proton at high Q 2 scales. 



I. INTRODUCTION 

The production of charm quarks in association with 
electroweak gauge bosons at hadron colliders can pro- 
vide important information on strange and charm quark 
PDFs, complementary to that obtained by tagging charm 
quarks in the final state in deep inelastic scattering ex- 
periments |l|. In particular, the Tevatron CDF and DO 
experiments @, Q have measured the cross section for 
charm quarks produced in association with W bosons, 
using muon tagging of the charm-quark jet. However the 
accuracy of these measurements is limited to ^30% by 
low statistics. The LHC is expected to provide a more 
precise measurement, and indeed the CMS collaboration 
has recently performed a similar study [J] of W^ + c(c) 
production, again using muons to tag the charm quark 
jet in the final state. 

At leading order (LO), the Feynman diagrams for 
charm production in association with a W boson are 
shown in Fig. [T] The dominant contribution comes from 
strange quark - gluon scattering, as the corresponding 
down-quark contribution is strongly Cabibbo suppressed. 
The cross section for W + c production (where 'c' denotes 
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FIG. 1: Feynman diagrams for W +c(c) production at lead- 
ing order. 

a tagged charm-quark jet) is measured in Ref. 2] (CDF), 
while Ref. [3] (DO) introduces the ratio of charm jets to 
all jets, which is expected to suffer less from both ex- 
perimental and theoretical uncertainties. CMS has per- 
formed a similar analysis [4] using the 2010 LHC data 
set. 

Also of interest is charm production in association with 
Z bosons. The leading-order process is simply eg — > Zc, 
and so this process can be used to extract informa- 
tion on the charm quark PDF. It is important to note 
that for both W + c and Z + c production at hadron 



colliders, the strange and charm quarks are probed at 
much higher Q 2 (~ 10 4 GcV 2 ) values than in the tra- 
ditional determinations from deep inelastic scattering, 
i.e. vs — > ji~ c{-^- /i + ) and ec — > ec(— > /i + ) where typ- 
ically Q 2 ~ 10°~ 2 GeV 2 ). Taken together, the mea- 
surements therefore also test DGLAP evolution for these 
quark flavours. 

In this letter we study VF + c-jet production in the con- 
text of the CMS analysis [4j , analysing the different quark 
contributions and comparing the predictions of various 
widely-used PDF sets. We also study the corresponding 
cross-section ratio for Z + c-jet production, which should 
be measurable with the 2011 LHC data set. 



II. CMS MEASUREMENT OF <r(W + c) 

The two relevant cross-section ratios introduced by 
CMS 4] are: 



R,. 



a{W + 



a(W-+c) 



and R c 
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a(W +jet) 



(1) 



The advantage of using ratios is that many of the the- 
oretical and experimental uncertainties cancel. In par- 
ticular, the ratios are fairly insensitive to higher-order 
perturbative QCD corrections. Note that Rf = 1 at 
the Tevatron. We calculate the cross sections at NLO 
pQCD using MCFM @, applying the CMS cuts [| to 
the final-state: p{, > 20 GeV, \rf>\ < 2.1, p l T > 25 GeV, 
|V | < 2.1, R = 0.5, R lj = 0.3. Five different NLO 
PDF sets are used: CT10 6], MSTW2008 [7], NNPDF2.1 
0, GJR08 |9j and ABKM09 [1Q], as implemented in 
LHAPDF 11]. The renormalisation and factorisation 
scales are set to Mw, although the cross-section ratios 
are rather insensitive to this choice. (We have also con- 
sidered dynamical scales of the form Q 2 = M^ + p^ w , 
but the differences for the cross-section ratios are similar 
in magnitude to the PDF uncertainties.) 

The results are summarised in Table U where we also 
include the ratio: 



CT (VF+ +jet) 
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o-(W~ +jet)' 



(2) 



Ratio 


R c 


Rc 


R ± 


CT10 


n Q5Q+0.009 
u.»uo_ 007 


124+ ' 021 


1 QQ+0.03 


MSTW2008NLO 


nq21 +0.022 
u.»zi_ 033 


11 fi+ooo 2 

U.110_g q 02 


1 34+0.01 
J -°*— o.ai 


NNPDF2.1NLO 


0.944±0.008 


0.104±0.005 


1.39±0.02 


GJR08 


0.933±0.003 


0.099±0.002 


1.37±0.02 


ABKM09 


0.933±0.002 


0.116±0.003 


1.39±0.01 



TABLE I: Comparison of results for the cross-section ratios 
defined in ((TJ and J2| at NLO using different PDF sets, in- 
cluding 68%cl (asymmetric, where available) PDF errors. 



We note that our values for Rf are systematically 
larger by ~3% than those quoted in the CMS study [4j, 
but we are unable to account for this difference. We 
have also checked that the values are stable with respect 
to NLO corrections, with the difference from the LO re- 
sults being within 1% for Rf and 3% for R c . In what 
follows we restrict our analysis to the three global PDF 
sets CT10, MSTW2008 and NNPDF2.1, since these span 
a relatively broad range of predictions for the cross sec- 
tion ratios. 

For reference, we note that the values of Rf and R c 
measured by CMS are: 



Rf = 0.92 ± 0.19 (stat.) ± 0.04 (syst.) 
R c = 0.143 ± 0.015 (stat.) ± 0.024 (syst.) 



(3) 
(4) 



Note that the experimental systematic errors are already 
of the same order as the differences between the predic- 
tions of the various PDF sets. 

If the strange quark contribution to W + c production 
were totally dominant, then any deviation of Rf from 1 
would imply an asymmetry between s and s. However 
even if s = s, the fact that d < d will automatically 
give Rf < 1 through the Cabibbo suppressed d— quark 
contribution. Schematically, at LO we expect 



R,, 



s + \V dc \ 2 d 
s+\V dc \ 2 d' 



(5) 



with Vd c =0.225. This leads to a suppression by a factor 
of 20 of the d— quark contribution. 

The contributions (in fb) from d and s quarks to the 
LO W + c cross sections are shown in Table [TT| and the 
corresponding percentages in Table ULTl The cross-section 
values correspond to the leptonic decay channel W —> ev. 
These are obtained using LO subprocesses (but calcu- 
lated with NLO PDFs) since the presence of additional 
NLO subprocesses does not permit a simple quark-flavour 
decomposition. 

The relative contributions of initial-state s and d 
quarks to the W + c cross-section ratios R^r and R c are 
illustrated in Figs. [2] and [3] For CT10 s = s and there- 
fore the fact that Rf < 1 is due entirely to the difference 
between d and d in the Cabibbo suppressed dg — > Wc 
subprocess. NNPDF2.1 does have an asymmetric strange 



Process 


CT10 


MSTW2008NLO 


NNPDF2.1NLO 


W + + c: s 


39934 


37133 


32980 


W + + c: d 


2666 


2854 


2880 


W~ + c: s 


39987 


38449 


33012 


W~ +c: d 


4969 


5178 


5180 



TABLE II: Cross section contributions (in fb) of d— quarks to 
a{W + c) ■ B{W ->• ev) for different NLO PDFs at LO. 



Process 


CT10 


MSTW2008NLO 


NNPDF2.1NLO 


W+ + c 
W~ +c 


6.3 
11.1 


7.1 
11.9 


8.0 
13.6 



TABLE III: Percentage contribution of d— quarks to a(W + c) 
for different NLO PDFs at LO. 



sea, s — s ^ 0, but the asymmetry turns out to be very 
small in the x, Q 2 region of interest for this process and 
therefore R^r ^ 1 is again determined mainly by the d, 
d asymmetry. Finally, for MSTW2008 the strange asym- 
metry is larger and therefore contributes significantly to 
Rf. Here the (LO) ratio of <r{W+ + c) to a(W~ + c) ob- 
tained by setting d— quark PDFs to zero is ~ 0.96, and 
this is decreased further by the asymmetry in d and d. 

The strange asymmetry s — s for Q — M\y, the relevant 
scale for this process, is shown in Fig. [H including the 
68%cl uncertainty band in the case of MSTW2008NLO 
and NNPDF2.1. The strange asymmetry in both of these 
sets is constrained by the CCFR and NuTeV dimuon 
vN and DN DIS data [II] in the global fit. These data 
slightly prefer an asymmetric strange sea in the x range 
0.03 — 0.3, although the CT10 symmetric choice of s = s 
is also consistent with the data within errors. In the 
MSTW2008NLO fit, the strange asymmetry is parame- 
terised 'minimally' as 



sv{x,Ql) 



s(x,Ql) -s(x,Ql) 



= A- 



„<5_- 



^(l-x^-il-x/xo), (6) 



with the overall constraint that L dx sy(x, Q 2 ,) = 0. It is 
this choice of parametrisation that drives the relatively 
large positive asymmetry in the range x ~ 0.01 — 0.1. 
There is no such strong parameterisation dependence in 
the NNPDF2.1 fit. A precise measurement of the ratio 
R^ , combined with an improved knowledge of the d, d 
difference (for example, from the rapidity dependence of 
the inclusive W — ¥ iv charge asymmetry), could there- 
fore provide important new information on the strange 
sea asymmetry at small x. 

The cross-section ratio R c can be used as a measure 
of the total strangeness of the proton, and to the ex- 
tent that these W+jet cross sections are dominated by 
qg scattering we can expect R c ~ w ~+% ■ F° r our three 
sets of NLO PDFs this ratio at scale Q = Mw is shown 
in Fig. \5\ We note that the ordering of the R c values 
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FIG. 2: Effect of initial-state s and d quarks on R^r using 
NLO PDFs and LO subprocesses. 
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FIG. 3: Effect of s and d quarks on R c using NLO PDFs (LO 
processes). 



for the different PDF sets, see Table |H agrees qualita- 
tively with the corresponding values of the quark ratio 
at x ~ 0.06, the average value of the incoming quark x 
for this collider energy and choice of cuts. We also see 
that the relative size of the PDF uncertainties in Table Q] 
aligns well with the PDF uncertainties of Fig. \5\ In par- 
ticular, the MSTW2008NLO strange-quark error band is 
much narrower than that of the other sets because of the 
implicit assumption in the MSTW global fit that all sea 
quarks have the same universal Qi{x, Qq) ~ x behaviour 
as x — > 0. In practice, the parameter 5 is determined 
quite precisely by the fit to the HERA small- x structure 
function data. 

Of the two ratios, Rf is less affected by the choice of 
selection cuts. Modifying the lepton and jet cuts, within 
experimentally reasonable ranges, changes R^r by only a 
few percent. R c is more sensitive to the selection cuts 
with the value changing by up to 0(30%). 

When in future more high-statistics data are available, 
the ratios Rf and R c can also be considered as distri- 
butions of kinematic observables, e.g. the W transverse 
momentum as shown at LO (using NLO PDFs) in Fig. [5] 
for Rf. In contrast to R^, which is related to the u/d 
ratio at high x and therefore increases with p^ , R^r is a 
decreasing function of p^ driven by the dominance of the 
valence d— quark at high x over the other parton distri- 



FIG. 4: Strange valence distribution for NLO PDFs at Q = 
Mw. For MSTW2008NLO and NNPDF2.1, the shaded bands 
correspond to the 68%cl PDF uncertainty. 
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FIG. 5: NLO PDF ratio of s + s to E(g + q) at Q = Mw- 

butions involved, see ([5]). This decrease of Rf is obtained 
for all PDF sets but the rapid drop for NNPDF is a re- 
sult of both d/d at large x and also the increasing value 
of sy at large x as shown in Fig. 01 Similar conclusions 
can be drawn by considering Rf as a function of the W 
rapidity as shown in Fig. [7] 

III. PREDICTIONS FOR a(Z + c) 

Even though the corresponding cross sections with the 
W boson replaced by a Z boson are significantly smaller, 
especially when account is taken of the difference in lep- 
tonic branching ratios, with a sufficiently large data sam- 
ple a similar analysis can be performed. 

We first consider the ratio 



R„ — 



a(Z + c) 



(7) 



<7(Z+jet)' 
where the c in the numerator here refers to either a c or a 
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Dependence of Bit on pYp 7 using NLO PDFs. 
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FIG. 7: Dependence of Rf on y w using NLO PDFs. 



cjet. Defining a similar set of experimental cuts: p J T > 20 
GeV, \rf\ < 2.1, v l T > 25 GeV , W\ < 2.1, R = 0.5, 
R 1 ' = 0.3 and 60 < m u < 120 GeV, gives the NLO ratio 
predictions shown in Table HVl now with the QCD scales 
set to M z . (Note that MCFM also includes the photon 
(7*) contribution to the Z+jet cross sections, consistent 
with the experimental analysis, but the number of photon 
events is strongly suppressed by the cut on the lepton pair 
invariant mass.) 



PDF set 


Re 


CT10 


0619+ 00032 


MSTW2008NLO 


0640+ 00014 

U.UO4:U_ 0016 


NNPDF2.1NLO 


0.0660±0.0013 


GJR08 


0.0611±0.0011 


ABKM09 


0.0605±0.0019 



TABLE IV: Comparison of flf NLO predictions for the dif- 
ferent PDF sets, with 68%cl PDF uncertainties. 
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FIG. 8: Charm quark fraction (c + e)/£(g + q) at Q — M z 
for NLO PDFs. 



In principle R^ provides direct information on the 
charm content of the proton, complementary to that ob- 
tained from DIS experiments via F£. 

We note that the differences between the predictions 
of different PDF sets are much smaller than for the 
strange quark distributions, presumably because in all 
these global fits the charm distributions arise pertur- 
batively from g — > cc splitting, with the small- a; gluon 
well determined from the HERA structure function data. 
This can be seen in Fig. [8] which compares the ratio of 
charm quarks to all quarks for the three PDF sets. When 
PDF errors are taken into account the three sets are very 
difficult to distinguish, and this is reflected in the simi- 
larity in the i?f predictions in Table HVl The use of R c as 
a PDF discriminator will therefore require a very precise 
measurement. 

We can also consider the (charm) charge asymmetry 
ratio: 



Rt(z) 



a(Z + c) 
a(Z + c)' 



(8) 



At first sight it might appear that Rf(Z) is automat- 
ically equal to 1, but this is only true if c = c in the 
initial state. This is indeed the case for all the PDF 
sets considered here, since as already noted the charm 
distributions are generated by charge symmetric g —)■ cc 
splitting. However this symmetry does not necess arily 
hold if we allow for an intrinsic charm component [13j . 
PDF studies incorporating intrinsic charm, see for exam- 
ple p], [lj| , suggest that it is probably a small effect com- 
pared to perturbatively generated charm, particularly at 
the small x values relevant to the LHC. In any case, any 
such intrinsic charm analysis is highly model dependent, 
and beyond the scope of the present study. We do note 
however that in the process qg — > Zq, which dominates 
Z+jet production over most of the kinematic range at 
the LHC, the 'g' is more likely to be positively charged 
than negatively charged, and the q— jet is therefore more 
likely to contain a jx + (say) than a /i~ . Hence there may 



well be a 'natural' charge asymmetry in the misidentified 
charm-jet background. 

Finally we can also compare the ratios: 



TjWZ 



a{Z + c) 
a(W + c) 



and R wz = 



(T{Z+jet) 
cr(W + jet)' 



(9) 



The NLO predictions for MSTW2008NLO are 0.045 and 
0.082 respectively, for the selection cuts described above 
for W and Z including leptonic decays. We can relate 
the ratio of these ratios to the ratio of s + s over c + c 
which can be read from Figs. [5] and [5] at the appropriate 
momentum fraction. Of course in practice it is more 
complicated, as multiple scales, momentum fractions and 
couplings are involved, but an estimate R wz / R^ z ~ 2 
can be extracted for this ratio and this matches well with 
the 0.082/0.045 ratio above. 



IV. CONCLUSIONS 



strange content of the proton, and in particular the asym- 
metry between s and s at small x and high Q 2 . We 
have also shown results for several differential distribu- 
tions that can in principle be used to provide additional 
information on the x dependence of the strange and anti- 
strange quark distributions. We also propose a measure- 
ment of the corresponding ratio for Z bosons, R z , which 
can be used as a measure of the charm content of the 
proton. 

Finally, we note that information on the strange quark 
content of the proton can also be obtained by comparing 
the total W and Z cross sections at the LHC, see for 
example the recent ATLAS study reported in [15|. This 
exploits the fact that the dependence on the strange PDF 
is linear for the W and quadratic for the Z: a(W) ~ 
cs+cs+..., o{Z) ~ ss+.... The method is complementary 
to the charm-jet tagging method discussed in the present 
study, and it will be interesting to compare the results 
from the two analyses when more precise data become 
available. 



We have investigated charm production in associa- 
tion with W and Z bosons at the LHC. We have pre- 
sented predictions relevant to the recent (and ongoing) 
CMS analysis for W bosons produced in association with 
a charm-quark jet. Use of a larger LHC data sample 
should lead to precise measurements of the ratios R c and 
R~ which can in turn provide useful information on the 
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